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ABSTRACT 

We study the compression driven magnetic reconnection in the relativistic 
electron-positron plasma. Making use of a 2.5D particle-in-cell code, we simulated 
compression of a magnetized plasma layer containing a current sheet within it. 
We found that the particle spectrum within the reconnecting sheet becomes non- 
thermal; it could be approximated by a power-law distribution with an index of 
-1 and an exponential cutoff. 

Subject headings: plasmas, magnetic fields, acceleration of particles 



1. Introduction 



Magnetic reconnection is one of the important processes for the particle acceleration and 
heating. The production of energetic particles has been extens ively studied in non-relativistic 
reconnection (ICargill 200 ll ; iDrake et al. 20051 ; iPritchett 20061 ) , mostly in the solar terrestrial 
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ray bursts ( Drenkhahn 200 j : Drenkhahn fc Spruit 20021 : Thompson 20061 ). Genuinely rela- 



tivistic reconnection occurs if the magnetic energy exceeds the plasma energy including the 
rest energy. Then dissip ation of the magnetic field inevitably leads to relativistic energies of 
the particles (kirk 2004T ). 



Relativistic reconnection typically occurs in electron-positron plasmas; in this case 
physics of the process is simplified considerably as compared with electron-ion plasmas where 
two significantly different scales are presented. Larrabee, Lovelace & Romanova (2003) cal- 
culated self-consistent equilibria of the relativistic electron-positron plasma in the vicinity 
of the magnetic X-point in the reconnecting current sheet. The particles are accelerated 
there by the electric field parallel to the X-line. It was found that the particle distribution 
function is described by a power law dn/dj oc 7 _1 with an exponential cutoff. Note that 
generally only a small fraction of magnetic energy is released in the vicinity of the X-point. 



-2- 



Reconnection in the X-point just allows the magnetic field lines to shrink thus releasing the 
magnetic energy. For example, in the classical Petschek model, most of the energy is released 
at slow shocks that stem from the X-point. Therefore in order to find overall particle distri- 
bution function, one have to study particle acceleration/heating in the whole reconnection 
region; the results will presumably depend on global geometry and boundary conditions. 
Particle-in-cell (PIC) simulations of spontaneous reconnection in an infinite, plane current 
sheet (Zenitani & Hoshino 2001, 2005, 2007, 2008; Jaroschek et al. 2004) confirm that a dc 
acceleration takes place around the X-point and that the particle energy spectrum in this 
region is roughly described by a power law with the power index of -1. However, the energy 
spectrum over the whole simulatio n domain is significantly steeper; it is approxim ated by 



the power law with an index of -3 (IJaroschek et al. 2004 ; IZenitani fc Hoshino 20071 ). 



This paper is aimed at the particle acceleration in the compression driven reconnection. 
Our study is motivated by observation (Lyubarsky 2003, 2005; Petri & Lyubarsky 2007) that 
reconnection could be an essential part of dissipation mechanism at the termination shock in 
the striped pulsar wind. Pulsars lose their rotational energy predominantly on generation of 
the Poynting dominated winds. Most of the energy is transferred in the equatorial belt where 
the sign of the magnetic field alternates with the pulsar period forming stripes of opposite 
magnetic polarity; such a structure is called a striped wind (e.g., review by Kirk et al. 2007). 
When the striped wind arrives at the termination shock, the plasma is strongly compressed; 
in the comoving frame, the compression ratio is very large, about the Lorentz factor of 
the upstream flow. Therefore the alternating magnetic fields are easily annihilated by the 
compression driyen reconnection. This c onjecture is supported by 1.5D PIC simulations 



(ILvubarskv 20051 : IPetri & Lvubarski 20071 ). 



In one-dimensional simulations, the magnetic energy is transformed into heat because 
all the particles gain energy with the same rate. In order to study particle acceleration, 
multidimensional simulations are necessary because one can expect a formation of non- 
thermal tails in particle spectrum if different particles gain different energy; this could happen 
if X-points are formed within the current sheet. As a step towards the multidimensional 
simulations of the shock in a striped wind, we performed 2.5D PIC simulations of driven 
magnetic annihilation within a plasma layer containing only two stripes of the opposite 
magnetic polarity. In our simulations, the plasma layer is compressed by an external force, 
which imitates the compression within the shock structure. We show that within such a 
structure, a non-thermal particle spectrum is formed, which can be roughly described by a 
power law with the slope -1 and an exponential cutoff. We analyze the particle motion within 
the compressed layer and show that the particle acceleration in the vicinity of X-points plays 
a crucial role even though only a small fraction of the energy is released there. The reason 
is that in the compressing medium, the larger the particle Larmor radius, the more energy 
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it gains therefore particles preaccelerated in the X-point gain more energy than the particles 
bypassed the X-point. Thus a large fraction of the total energy is eventually transferred to 
particles that passed the X-point. 

The paper is organized as follows. In Section 2, simulation parameters and methods are 
introduced. In Section 3, the results are shown and discussed. Section 4 gives the conclusions 
of this work. 



2. Simulation setup 

We used 2.5D (2D3V - 2 spatial and 3 velocity components) fully relativistic electro- 
magnetic particle-in-cell code. The evolution of electric and magnetic fields is governed by 
Maxwell equations for components of the field vectors E x , B y and B z . In two-dimensional 
configuration with B x = 0, the components E y and E z of the electric field are decoupled 
from the remaining field components and could arise only due to charge fluctuations (cf. 
simulations by Karlicky (2007) and Zenitani & Hoshino (2008) where the charge separation 
arises at B x ^ 0). In electron-positron plasma with the same distribution functions for both 
species, such fluctuations are very low and therefore we take E y = E z = throughout the 
simulations. Such a suppression of the electrostatic fluctuations permits the use of relatively 
small number of particles. 

The fields are u pdated usi ng the leap-frog scheme. The staggered grid mesh system, 



known as Yee lattice (lYee 19661 ). ensures that the change of the magnetic flux through a cell 
surface equals the negative circulation of the electric field around that surface and the change 
of the electric flux through a cell surface equals the circulation of the magnetic field around 
that surface minus the current through it. Then the divergency-free condition is maintained 
to the machine accuracy. Here the electric and the magnetic fields are in a symmetry form 
except subtracting the charge flux J in Ampere equation. The curr ent density is ca lculated 



and subtracted after the particles are moved later in the program (IBuneman 19931 ). Before 
and after moving (or pushing) the particles, the magnetic field is updated in two half steps 
so that it is available at the same time as the electric field for the particle update. The 
particles positions and velocities are advanced by Newton-Lorentz 3D equat ions of motion, 



which were solved by the Buneman-Boris method (IBirdsall fc Langdon 19851 ). 



Assuming the reconnecting current sheet to lie in the y-z plane, FigJTh. illustrates the 
configuration of the simulated domain (rectangular grid shown by dashed line). The system 
lengths in the y- and z-directions are L y and L z respectively. The domain is divided into 
3000(L<,) x 400(L^) cells grid or normalized to electron skin depth 384(LyUJ p /c) x 51(L z lu p /c), 
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where the plasma frequency uj p is calculated for the particle density in the lobe. Fig (lb 
demonstrates the initial profiles of density n(y), temperature T(y) and velocity Vd r (y) near 
the current sheet (y = 1400 — 1600). Outside the central part of the simulation box, the 
density remains constant up to the box margins < y < 10 and 2990 < 3000. 




Fig. 1. — (a) Schematics of initial configuration, fields and charge fluxes, (b) Initial profiles 
of density n(y), temperature T(y) and velocity Vd r (y) near the current sheet. Outside the 
central part of the simulation box, the density remains constant up to the box margins 
< y < 10 and 2990 < 3000, where the plasma is absent (n = 0) 

To avoid escaping in the y-z plane, the particles are seeded away from boundaries at a 
distance larger than particles gyroradius. Ten mesh points are set to keep the margin on each 
side of the boundaries (shown gray), i.e., particles are loaded in the region of 10 < y < 2990. 
At the boundaries, two external currents normal to the simulation plane are set (as illustrated 
in FigfT]by solid lines at gray margins) in the same directions. The magnetic field is directed 
along z\ in the region between the currents it is set initially as 



B z (y) = fiotanh 



L y -2y 



2A ' w 
where A(= 25) is the sheet halfwidth. The initial magnitude of the external currents is 

J« = ^ (2) 
Q 

so that the magnetic field vanishes in the region outside the currents on both sides. 

The magnetic field ([1]) implies the current in the x-direction; the current density is 
determined by Ampere's law 

VxB = ^. (3) 
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This current is generated by motion of electrons and positrons in opposite directions with 
the velocity 

c dB z (y) 

v dr{y) = , . -, — • (4) 
2qn(y) dy 

The initial spatial distribution of particles is non-uniform along ^-direction 

/ s °ft"max ^0 i r \ 

" iV> = "° + co 3h ((V2- y )/A) ^ (5) 

In the main simulations presented here, the total number of particles is 5 • 10 6 , giving no = 4 
particles per cell in the lobe and n max = 10 in the sheet. In order to check reliability of our 
simulations, we have also run similar models with smaller box 300 (L y ) x 40 (L z ) but with 
different number density of particles, where n was multiplied by factor a = 1 -f- 10. In order 
to keep the system in mechanical equilibrium (jHJ), the magnetic field B was multiplied by */a, 
while other parameters were kept unchanged. We found that particles spatial distribution 
and energy spectrum during the compression were independent of the number density n. 

The particles are initialized by the Maxwell distribution with a temperature depending 
on y such that the system is in mechanical equilibrium, i.e. 

B 2 

nT H = const. (6) 

87T 

Then 

T(y) = B o~ B M 2 + Tq (7) 



167rn(y)yl - (v dr (y)/c) 



where To is a minimal temperature value, To = 0.1m e c 2 . The chosen initial state is close to 
the true kinetic equilibrium because the initial thickness of the sheet is significantly larger 
than the initial particle Larmor radius, 2A/r^ w 29, where = , ymvc/qB w 1.75 (the 
thermal velocity v is calculated using the initial temperature at the center of the sheet, 
which is T = 1.22mc 2 ). 

The imposed boundary conditions along z-axis are periodic. Along y-axis, we impose 
Landman's radiation-absorbing condition, which requir es that the field s should be able to 



radiate away into space and should not be reflected (ILindman 19751 ). The time step is 
chosen such that u c At — 0.2 where uj c = qBo/mc. The plasma in the lobe is magnetically 
dominated; uj p /uj c is initially 0.32 (no is used to calculate u p ) and decreases in the course of 
compression. 
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3. 



Simulation results 



In simulations, the reconnection is driven by compression. We increase the magnitude 
of boundary currents on both sides so that the magnetic field grows and the plasma is 
compressed towards the midplane of the simulated domain. Note that in the course of com- 
pression, the plasma layer becomes narrower so that the pressure balance is restored faster 
therefore in order to save the computer time we can take the compression rate increasing 
with time. We choose the boundary currents quadratically growing with time: 



where v = 1.3. We have also run the simulation with the compression rate v = 0.75 (twice 
slower) however no significant difference in the final spatial distribution and energy spectrum 
was observed. 

As a result of external compression, the magnetic field lines reconnection occurs in the 
midplane where the field reverses sign. Fig J2] demonstrates evolution of spatial distribution of 
particles. The simulated domain is initially separated into several stripes, which are located 
at different distances from the midplane. The particles located initially in each stripe are 
marked with a corresponding color which allows tracing the motion of different pieces of 
plasma in the course of compression. FigEti shows initial distribution of particles, where due 
to the reflection symmetry, the stripes at the same distances from the midplane are colored 
identically on both sides. The width of the stripes is chosen the same except of the central, 
more dense, stripe (cyan) which is twice as large as the sheet width. As the boundary 
currents grow, the plasma is compressed and each stripe shifts towards the midplane i.e. to 
y = 1500, (see FigfSb). In the lobe region, the particles energy is low, therefore the Larmor 
radii are small and particles are restricted to move with the magnetic field lines towards 
the current sheet. Therefore each stripe moves as a whole until it reaches the current sheet. 
Within the current sheet, the magnetic field lines are reconnected and as a result particles 
assemble on the magnetic island, where each stripe forms a corresponding ring, as one can 
see in Figj2U. 

Evolution of the particle energy spectrum in each stripe is presented at Figj3j Initially, 
as shown on Figj3h., the energy of particles in the lobe is low and only in the current sheet, 
the higher energy particles are present. When compressed, the stripes begin to form ring- 
regions in the sheet. When a stripe assembles in a ring, the particle energy spectrum extends 
to higher energies FigJHfc shows the spectrum of the magenta stripe, which begins to form a 
ring at t — 7000. Figj3H shows that by the end of simulations, the hard tail in the energy 
spectrum is formed predominantly by particles entered the current sheet in the course of 



J(t) = Jo 1 + v—t 





y 

Fig. 2. — Evolution of the particle spatial distribution initially separated into colored stripes 
along y-direction at (a) t = 0, (b) t = 5000, (c) t = 7000 and (d) t = 9990. 

the reconnection process. The fraction of particles residing in the sheet from the beginning 
(cyan curve) is small therefore the final spectrum is independent of the initial configuration 
of the system. 

In the course of compression and magnetic reconnection, all the particles gain energy. 
However, taking two identical particles at the same distance from the sheet center one can 
find a significant difference (by two orders of magnitude) in their final energies. This wide 
range of final energies demands further explanation. We take a number of particles with high 
and low final energies initially located at the same distance of 300 from the midplane (i.e. 
y = 1200 and y = 1800) and trace them during the compression process. FigJH demonstrates 
the energy evolution of high-final (solid curves) and low-final (dashed curves) energy test 
particles. As shown, the high-final energy particles begin acceleration at approximately the 
same time (t ~ 5000) when the chosen stripe enters the current sheet. The particles from 
the same stripe that were not accelerated at this time remain cold. In Figj5j five typical 
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Fig. 3. — Evolution of the particle energy spectrum for each of the colored stripes at (a) 
t = 0, (b) t = 5000, (c) t = 7000 and (d) t = 9990. 

trajectories of high-final energy particles are shown during the short time-interval (At = 100) 
when the particle energy just begins to grow. The simultaneous magnetic field lines are 
shown by thin lines. One sees that, during the onset of acceleration, all these particles are 
located in the vicinity of the X-point. Afterwards the particles leave the X-point however 
their energy continues to grow (see FigfjJ, whereas the energy of particles that bypassed the 
X-type region remains low. The reason is that in the course of compression, the particles 
gain energy proportionally to the available one. Therefore the particles preaccelerated near 
the X-point continue to gain energy while those that bypassed the X-point remain cold. 

A closer look at the particle trajectory shows that the steepest jumps in energy appear 
at certain locations within the X-point region. The particle energy grows sharply when 
the electric field becomes equal or even exceeds the magnetic field. In order to trace the 
acceleration sites, the quantity 

B 2 — E 2 
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Fig. 4. — Energy of test particles. Solid - high-final energy particles. Dashed - low-final 
energy particles. 

was constructed and shown in FiglHk as a color background; here x ~~ * — 1 corresponds to 
E ^> B and \ corresponds to E <C B. The locations of the test particles are also 

presented for high-final (colored dots) and low-final (yellow crosses) energy particles. FigJ6b 
shows the energy evolution of the traced particles; the vertical dashed line marks the instant 
(t = 8260) of the snapshot at Figj6^. One can see that when a test particle (the green and 
blue dots in FigJHH) enters the region of low %-value, the energy (the green and blue curves 
in the bottom panel) of the particle explosively grows. 

The black and cyan particles are already outside the X-point region; they were acceler- 
ated in this region earlier and now gain energy by compression. One sees in Figj6b that the 
energy of the cyan particle varies in oscillatory manner; at the later stage, the green particle 
exhibits the same oscillatory behavior. In order to explain such oscillations, the electric and 
magnetic fields at the same time t = 8260 as Figj6^ are shown in Fig. [7^ and b, respectively. 
One sees a small magnetic island in the region where the cyan particle is located in Figj6^i. 
The island rapidly moves away from the X-point region as shown in Fig. [7b by arrows and 
eventually merges into the large island at the upper half of the current sheet. Such islands 
are formed occasionally when the magnetic field lines are teared in the X-point. The electric 
field within such a rapidly moving island reverses sign according to E = v x B, being posi- 
tive in front and negative in back of the island motion, as one can see in Fig. [TJi. Note that 
within the island, the electric field remains smaller than the magnetic field (even though the 
magnitude of the electric field is significantly larger within the island than near the X-point) 
therefore the particles are not accelerated there. Low energy particles just drift together with 
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Fig. 5. — Magnetic field lines (gray) and five typical trajectories (black) of high-final energy 
particles during the onset of acceleration at t ~ 5000 (see Fig. H]). All the trajectories (over 
the interval At = 100) are within the X-type region. 



the island because their Larmor radius is less than the size of the island however particles 
with the Larmor radius comparable with the size of the island experience oscillations. 

When a small magnetic island appears, an additional X-point arises; in Fig. [7b it is at 
the top of the box (recall that the box is periodic in the ^-direction). An additional low-x 
region arises there (at the top of Fig. Ek) so that an additional acceleration cite appears. 
In Fig. [6^, the red particle moves down so that it is going to enter this region from above; 
Fig. [6b shows that this particle is indeed accelerated soon. 

FigJH]shows evolution of the overall energy spectrum. One sees that the range of particles 
energies gradually expands and the high-energy tail appears as plasma is compressed. The 
resulted energy spectrum is well approximated by the expression proposed by Larrabee et 
al. (2003): 

— oc - exp ( — - I . 10 

«7 7 V 7o/ 

They showed that this formula fits the energy spectrum in the vicinity of the X-point. In our 
simulations, the overall particle distribution is well fitted at high energies by Eq. ffTOj) . The 
reason is that the high energy tail is formed only by particles preaccelerated in the X-point. 
Outside the X-point, they gain energy predominantly by compression. In this case the final 
particle energy is proportional to the initial energy so that the particle energy spectrum 
formed near the X-point just shifts towards larger energies, the spectral shape remaining 
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Fig. 6. — (a) Low (yellow crosses) and high (color) energy particles position over %-value 
background at t = 8260, (b) Energy evolution of the traced particles. 

unaltered. 

Finally, in order to show that the above magnetic reconnection is driven solely by the 
compression, the same simulation was performed for the uncompressed case with constant 
boundary currents on both sides. The timescale of the simulation was the same as in the 
compressed case, t = 9990, which corresponds to ~ 3.5 light travel times from the boundary 
to the center of the box (note that in the lobe the plasma is magnetically dominated therefore 
Alfven's velocity is close to the speed of light). Fig. [9] shows the particles final (at t = 
9990) spatial distribution, corresponding energy spectra and magnetic field. As one can see 
compared to Figs. [2^ and, [3^, practically nothing has changed as compared with the initial 
state; the particles remain at their initial locations along the y-direction (see Fig. (9H) and 
the energy spectra in each stripe have no considerable changes since the initial state (see 
Fig. [9b)- The small magnetic field islands observed in the current sheet midplane (Fig. [9b) 
develop already at the early stages of simulation whereupon nothing changes. The magnetic 
energy in the box also remains constant, to within the error of computation, throughout the 
simulation. These results clearly show that the initial configuration is stable, which means 
that the magnetic reconnection and the resulted particles acceleration discussed in present 
paper occur exclusively due to external compression. 
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Fig. 7. — Electric (a) and magnetic (b) fields in the current sheet at t — 9820 

4. Conclusions 

We performed 2.5D PIC simulations of the particle acceleration in the course of the com- 
pressionally driven reconnection in relativistic electron-positron plasma. Since the plasma is 
strongly compressed at the front of a relativistic shock, we consider our simulations as a sim- 
ple model for the particle acceleration at the termination shock in the striped pulsar wind. 
We found that the X-points play crucial role even though a relatively small fraction of the 
total energy is released there. Particles preaccelerated in the vicinity of the X-point take a 
significant fraction of the released energy because in the compressing medium, particles with 
larger Larmor radii gain more energy. This result could be of more general nature because 
for particles with larger Larmor radii, the frozen-in condition is violated easier therefore 
one can expect that in various situations, particles preaccelerated in the vicinity of X-points 
could take a good fraction of the dissipated magnetic energy. 

In our simulations, the fraction of the particles initially confined in the current sheet is 
small so that in the final state, the energy spectrum of the particles in the sheet is independent 
of the initial condition. This spectrum is found to be approximately dn/d^ oc 7 _1 for 
7 < 7o, while is rapidly decreasing for 7 > 7o- The maximal energy 70 grows in the course of 
reconnection. The electrons with the spectrum dn/d r y oc 7" 1 emit synchrotron radiation with 
the flat spectrum. Therefore our results support the idea that flat radio spectra of plerions 
may be attributed to the particle acceleration at the termination shock in the striped pulsar 
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Fig. 8. — Evolution of the total particle energy spectrum. The dotted curves represent the 
fit ( DUD - The spectrum curves are presented for the time t = 4000, 6000, 8000, 10000 and fit 
7o = 4, 16, 25, 35, correspondingly 

wind (Lyubarsky 2003). 

Note that the very occurrence of X-points was questioned by Zenitani & Hoshino (2005, 
2007) who showed that in the electron-positron plasma, the tearing instability, which is 
thought to be responsible for formation of X-points, develops slower than the drift kink 
instability, which just shifts straight magnetic field line tubes making the current sheet 
folded. Their PIC simulations show that the field dissipation due to the drift-kink instability 
does not result in non-thermal particle acceleration, the plasma is just being heated. The 
reason is that in this case, the field line tubes remain straight so that all particles gain energy 
with the same rate. However, the drift kink instability is suppressed in the presence of the 
current-aligned magnetic field (the so called "guide field"); then the magnetic reconnection 
occurs in X-points and produces a lot of nonthermal particles (Zenitani & Hoshino 2008). 
Therefore one can believe that the mechanism identified in this work remains relevant at 
quite general conditions. 

This work was supported by the German-Israeli foundation for scientific research and 
development under the grant 1-804-218.7/2003. 




-14- 




Fig. 9. — The uncompressed configuration at t — 9990: (a) particle spatial distribution, (b) 
energy spectrum for each of the colored stripes and (c) magnetic field lines in the current 
sheet. 
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